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ABSTRACT 

We present the results of a study of the globular cluster systems of 6 massive spiral galaxies, 
originally cataloged as low surface brightness galaxies but here shown to span a wide range of central 
surface brightness values, including two intermediate to low surface brightness galaxies. We used the 
Advanced Camera for Surveys on board HST to obtain photometry in the F475W and F775W bands 
and select sources with photometric and morphological properties consistent with those of globular 
clusters. A total of 206 candidates were identified in our target galaxies. From a direct comparison 
with the Galactic globular cluster system we derive specific frequency values for each galaxy that 
are in the expected range for late- type galaxies. We show that the globular cluster candidates in 
all galaxies have properties consistent with globular cluster systems of previously studied galaxies 
in terms of luminosity, sizes and color. We establish the presence of globular clusters in the two 
intermediate to low surface brightness galaxies in our sample and show that their properties do not 
have any significant deviation from the behavior observed in the other sample galaxies. Our results 
are broadly consistent with a scenario in which low surface brightness galaxies follow roughly the same 
evolutionary history as normal (i.e. high surface) brightness galaxies except at a much lower rate, but 
require the presence of an initial period of star formation intense enough to allow the formation of 
massive star clusters. 

Subject headings: galaxies: spiral — galaxies: star clusters — globular clusters: general 



1. INTRODUCTION 

Our understanding of the distribution of the central 
surface brightness of late-type galaxies has evolved signif- 
icantly in the last decades. Late-type galaxies were orig- 
inally thought to show a limiting central surface bright- 
ness value in the B band of /Lto=21.65 mag/arcsec 2 (the 
so-called Freeman value; Freeman 1970), but Disney's 
(1976) realization of the severe selection effects that the 
sky brightness introduces in the observations led to the 
current knowledge of the existence of significant num- 
bers of disk galaxies with much lower values of cen- 
tral surface brightness. McGaugh, Schombert & Bothun 
(1995) showed, using a selection-effects-corrected surface 
brightness distribution, that the number of galaxies faint- 
wards of the Freeman value remains approximately con- 
stant. This implies that the spatial density of the ob- 
jects known as "low surface brightness galaxies" (no > 23 
mag/arsec 2 ) is ~10 5 times higher than what is expected 
based on a distribution based on Freeman's value. 

Low surface brightness (LSB) galaxies encompass 
many of the extremes in galaxy properties. The vast 
majority are late-type disk galaxies with high gas mass 
fractions (Schombert et al. 1992), the surface density of 
this gas being well below the critical threshold density 
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for star formation (van dcr Hulst ct al. 1993). Therefore, 
as it has been confirmed from Ha studies (Schombert 
1990), LSB galaxies have very low star formation rates. 
They also appear to be more isolated than normal high 
surface brightness (HSB) galaxies (Bothun et al. 1993, 
Mo et al. 1994, Rosenbaum & Bomans 2004) and typi- 
cally show low metal content (de Blok & van der Hulst 
1998). All this observational evidence seems to indicate 
that LSB galaxies are relatively unevolved and quiescent 
objects, and may therefore provide an insight into the 
evolution of galaxies in unperturbed environments (de 
Blok et al. 1995). 

Several evolutionary scenarios have been proposed to 
account for the properties of LSB galaxies. Based on 
their observed blue colors, de Blok et al. (1995) con- 
clude that LSB galaxies cannot be the faded remnants 
of normal disk galaxies. More recently, van den Hoek 
et al. (2000) have shown that LSB galaxies roughly fol- 
low the same evolutionary history as HSB galaxies, ex- 
cept at a much lower rate. Using near-infrared observa- 
tions of a sample of 88 LSB galaxies, Galaz et al. (2002) 
show that a high fraction of LSBs have a well-developed 
old stellar population and that older LSBs are more fre- 
quent than optical data suggests. These observations 
are consistent with the scenario of Dalcanton, Spergel 
& Summers (1997), in which LSB galaxies are formed 
within a hierarchical formation scenario from low mass 
and/or high angular momentum proto galaxies which 
naturally form low baryonic surface density disks (see 
also Jimenez et al. 1998). On the other hand, Zackris- 
son, Bergvall & Ostlin (2005) argue using observations 
of extremely blue LSB galaxies that their properties are 
inconsistent with constant star formation rates over cos- 
mological timescales and that current observations can- 
not rule out the alternative possibility that these objects 
formed as recently as 1-2 Gyr ago. Therefore, some LSB 
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galaxies could have formed very recently. 

Globular cluster (GC) systems around galaxies are use- 
ful tracers to constrain the formation and evolution of 
their hosts (e.g., West et al. 2004). They trace all the ma- 
jor epochs of star formation in a galaxy and their lumi- 
nosities and colors contain information about the chem- 
ical enrichment environment and epoch in which they 
formed. Hence, the identification and study of the prop- 
erties of GCs in LSB galaxies will allow to obtain a phys- 
ical understanding regarding the differences and similar- 
ities between the formation and evolution processes of 
LSB and normal galaxies by comparing the properties of 
their GC systems. 

Here we present the results of a study of GC systems 
in a group of 6 galaxies, originally selected as LSB that 
turned out to span a range in central surface brightness, 
from normal to low surface brightness ones. This is the 
first study that aims to identify such objects in massive 
LSB galaxies. A previous study by Sharina et al. (2005) 
identified a number of cluster candidates in LSB galaxies 
using WFPC2 images, but their sample was constrained 
to dwarf galaxies. Our aim in this work is to probe for 
the presence of globular clusters in massive LSB galaxies 
and to compare their general properties to those observed 
in normal HSB galaxies. In particular, we aim to study 
the properties of objects consistent with being old GCs 
and which thus probe the early-stages in the formation 
of LSB galaxies. The mere presence of old GCs in these 
galaxies directly implies that there was an early period 
of star formation intense enough to produce massive star 
clusters. 

This paper is organized as follows. The observations 
and data reduction procedures are described in A 
study of the light profiles of the target galaxies aimed to 
obtain an improved classification of the galaxies in terms 
of their central surface brightness is presented in !j3l while 
^describes the selection of bona-fide GC candidates and 
a set of control-fields for each galaxy. In |JS] we use our 
observations to estimate the total number of GCs in our 
target galaxies and we discuss the properties of the de- 
tected GC systems in $6] This work is summarized and 
its implications are presented in EjTj 

2. OBSERVATIONS AND DATA REDUCTION 
PROCEDURES 

2.1. Selection of the Galaxy Sample 

Our sample of galaxies was selected from the LSB 
galaxy catalog of Bothun et al. (1985). Their work, based 
on the Uppsala General Catalog (UGC, Nilson 1973), 
classifies a galaxy as LSB if its average surface bright- 
ness /z pg = m pg + 5 log(Z?) + 8.63 (where D is the galaxy 
diameter in arcmin and m pg is the photographic magni- 
tude, both from the UGC) satisfies /i pg > 25 mag/arcsec 2 
(Bothun et al. 1985; O'Neil et al. 2004). This defini- 
tion is typically equivalent to the more commonly used 
classification of a galaxy as an LSB if the central sur- 
face brightness of the disk in the B-band, /xo,di s k, satis- 
fies /Uo jC iisk — 23 mag/arcsec 2 . The uncertainties in go- 
ing from the equation which describes a mean surface 
brightness to /io.disk are large though, and many galaxies 
satisfying the above criteria on // pg can have /io,disk w 22 
mag/arcsec 2 (O'Neil et al. 2004). We will come back 
to this issue in §3 below in order to clarify the general 



properties of our sample. 

We restricted this study to galaxies from Bothun et 
al. (1985) selected by their rotational velocity (mea- 
sured through the 21-cm velocity width W20) to have 
baryonic masses comparable to that of the Milky Way 
(Mmw ~ 7x 10 11 M Q ). As shown below, this implies 
that these systems should have GC systems rich enough 
to be detected in our targets and provide information on 
their ensemble characteristics to contrast with the wealth 
of information on the properties of GC systems in HSB 
galaxies. The sample was further selected to span a range 
of ~ 2 magnitudes in surface brightness and we applied 
a redshift cut-off (cz < 3000 km sec -1 or (m — M) < 33 
mag) to ensure that a large fraction of the GCs would 
be resolved using the ACS instrument on-board HST. 
Among the galaxies that satisfied these criteria, we se- 
lected 6 galaxies as a minimum number spanning the typ- 
ical ranges of surface brightness, H I mass, and Mm/ Lb 
ratios found in LSB galaxies. 

It is well known that the number of GCs in HSB galax- 
ies scales with the luminosity of the galaxy (e.g., Harris 
2001). This scaling is usually quantified through the spe- 
cific frequency Sm (Harris & van den Bergh 1981), which 
measures the number of GCs per unit U-band luminosity 
and is defined as S N = N GC 10°- 4 ( M ^+ 15 ), where M v is 
the total visual magnitude of the galaxy. The specific 
frequency is observed to depend on Hubble type, with 
early-type systems having Sn ~ 4 and late-type sys- 
tems Sn ~ 1, although the scatter around these values 
is significant. As no census of GCs in LSB galaxies ex- 
ists, we conservatively assumed a low specific frequency 
of Sn = 1 in predicting the expected number of GCs in 
LSB galaxies. With this assumption, we expected be- 
tween 50 and 150 GCs in our target galaxies. 

A list of our target galaxies and their basic properties 
is presented in Table [1] and a short discussion on the main 
features of each galaxy is included as an appendix. The 
distance to the galaxies was derived from the recession 
velocity corrected by Virgo infall using the local velocity 
field model given in Mould et al. (2000) with the Hub- 
ble constant fixed to Hq — 73 ± 5 km/s/Mpc, and were 
obtained from the Nasa Extragalactic Database (NED). 
Absolute magnitudes (M^J were derived with the most 
recent data compiled in NED and using the Schlegel et 
al. (1998) reddening maps. 

2.2. Observations 

Each galaxy was observed for a single HST orbit as 
part of program GO-10550 (PI: M. Kissler-Patig). The 
Wide Field Channel of the Advanced Camera for Sur- 
veys (ACS; Ford et al. 1998) was used to obtain five ~ 
400s exposures: three in the F475W band (« Sloan g) 
and two in the F775W band (~ Sloan i), which result 
in total exposure times of ^1200 sec and ~800 sec in 
the F475W and F775W bands respectively. In order to 
remove chip defects and bad pixels a line dither with a 
spacing of 0'.'146 in the g-band and 0'.'17 in the i-band 
was performed in the identical exposures of each filter. 
A log of the observations is given in Table [2] 

The choice of filters was dictated by the long baseline, 
which results in good sensitivity to metallicity and age 
of stellar populations, and by having at least one filter in 
common with the ACS Virgo and Fornax Cluster Surveys 
(Cote et al. 2004; Jordan et al. 2007a). These surveys 
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TABLE 1 
Properties of the galaxies." 



UGC00477 UGC03459 UGC03587 UGC06138 UGC11131 UGC11651 MW 



a( J2000) 


00:46:13 


5(J2000) 


+19:29:22 


I 


121.23 


b 


-43.36 


z 


0.008836 


v r [km/s] 


2698 


(m - M) 


32.84 


d [Mpc] 


36.98 


M B 


-19.14 


A v 


0.120 


E(B-V) 


0.036 


Mm ps [mag/arcsec 2 ] 


26.30 


Mo,disc[mag/arcsec 2 ] 


21.86 


W20 [km/s] 


246 


i [°] 


13 


^25 ['] 


3.5 


logfMm/iW©) 


9.85 


Iog(M dyn /M Q ) 


12.45 


Morphological Type 


Sdm 



06:23:59 06:53:54 11:04:39 

+04:42:39 +19:17:59 +27:43:26 

205.59 195.87 205.28 

-3.92 9.19 66.34 

0.009587 0.004226 0.008586 

2753 1227 2688 

32.88 31.13 32.83 

37.67 16.83 36.81 

-19.10 -18.33 

1.860 0.306 0.115 

0.561 0.092 0.035 

25.80 25.50 25.90 

21.43 22.63 23.74 

311 220 224 

53 22 31 

1.2 3 2 

9.05 9.30 9.43 

12.05 11.74 12.18 

Scd S? Sm 



1 8-1 0-99 


90-^7-1 ^ 




j-01 •s t v ? « 


+25-58-1 3 




29.64 


71.13 




9.85 


-12.58 




006014 


005087 




1958 


1730 




32.14 


31.87 




26.79 


23.66 




-19.58 


-19.06 


-20.5 


1.610 


0.879 




0.486 


0.265 




24.00 


26.50 




21.62 


21.75 




206 


276 




26 


15 




1.4 


3.5 




8.95 


9.24 


~9.7 


11.71 


11.97 


11.77 


Scd 


Sdm 


Sb-Sc 



a Positions, inclinations (i), 21cm velocity width (W20) and optical diameter at the 25mag/acrscc 2 isophote (D25) were extracted from Karen 
O'Ncil's Catalog of Massive LSB Galaxies (http://www.gb. nrao.edu/^koneil/biglsbgs/) and were used to derive the masses (Mjfj, WLdyn)- The 
photographic surface brightness /x m «« was computed from Uppsala General Catalogue (Nilson, 1973). The central disk surface brightness in the B 
band /io,disc was derived using the data presented in this work as described in the text. All the other parameters came either from NED or from 
HypcrLeda databases. Distances and distance moduli are corrected by Virgo infall. 



TABLE 2 
Observing log for GO-10550 



Galaxy 


Exposure Time 
F475W 


Exposure Time 
F775W 


UGC 00477 


3 x 404 


2 x 405 


UGC 03459 


3 x 400 


2 x 400 


UGC 03587 


3 x 404 


2 x 405 


UGC 06138 


3 x 406 


2 x 407 


UGC 11131 


3 x 400 


2 x 400 


UGC 11651 


3 x 406 


2 x 407 



observed 143 early-type galaxies and their GC systems 
in the Virgo and Fornax clusters using F475W (rj Sloan 
g) and F850LP (w Sloan z) and therefore matching one 
of their filters can prove useful to contrast the properties 
of our GCs with those of early-type galaxies. We chose 
F775W instead of F850LP due to its higher throughput. 
Note that here and throughout, we use g as shorthand 
to refer to the F475W filter, and i denotes F775W. 

2.3. Data Reduction and Object Detection 

The images were combined and cosmic-ray cleaned 
using the task "multidrizzle" in Pyraf (Koekemoer et 
al. 2002) without performing sky subtraction and by 
adopting a "lanczos3" kernel function for the final image 
combination. In order to register the images the header 
information was used, after checking using foreground 
stars in our frames that the actual shifts between the 
images were consistent with the commanded ones. The 
output of the image combination is an exposure time 
weighted image of 4218 x 4243 pixel 2 in each filter for 
each galaxy. Figure [1] shows the co-added F475W image 
for the six galaxies studied in this work. 

Object detection was carried out with SExtractor 
(Bcrtin & Arnouts 1996). In order to aid in the object 
detection and background determination a weight im- 
age was created using the WHT extension of the "mul- 



tidrizzle" final image, following the procedure outlined 
by Jordan et al. (2004). The result of this process, a 
map of the rms in the image, was used when running 
SExtractor by setting the "WEIGHT_TYPE" parame- 
ter to "MAP_RMS" . Object detection was performed in 
each image independently requiring a minimum of 5 con- 
nected pixels above a threshold of 2a. The background 
was estimated using a mesh size of 40 x 40 pixel 2 with 
no background filtering. The resulting object catalogs 
in each filter were matched using a matching radius of 
O'.'l and only sources that were matched were retained 
for further analysis. 

Aperture photometry was performed using an aperture 
radius of 4 pixels, with the background estimated locally 
in a rectangular annulus with a width of 30 pixels. The 
apparent AB magnitudes of the objects were calculated 
from the instrumental magnitude using the photomet- 
ric zeropoints and aperture corrections given by Sirianni 
et al. (2005) 7 . The dust maps of Schlegel et al. (1998) 
and the extinction ratios listed in Sirianni et al. (2005) 
were used to correct the observed magnitudes by fore- 
ground Galactic absorption. No corrections were applied 
to account for possible reddening internal to the target 
galaxies. 

7 The distances to our target galaxies imply that many of the 
GCs should be marginally resolved given their typical half-light 
radii ~ 3 pc (e.g., Jordan et al. 2005). This implies that the 
aperture corrections for point sources will in general underestimate 
the correction needed for GCs. We have checked the magnitude of 
the expected difference by creating simulated GCs as point-spread 
function convolved King (1966) models with half-light radius = 
3 pc and concentration c = 1.5, i.e. a model appropriate for a 
typical GC. We find that the average expected corrections vary by 
< 0.02 mag in g and < 0.05 mag in i. Given the modest size of 
the expected corrections for an average object we have chosen to 
adopt the corrections listed in Sirianni et al. (2005) for all target 
galaxies. Although such offsets do not affect our conclusions in 
any significant way, it should be kept in mind that magnitudes 
and colors can have systematic offsets < 0.05 mag due to the fact 
that GCs are marginally resolved. 




Fig. 1. — F475W band images of the six observed galaxies. From left to right, upper panel: ugc00477, ugc03459, ugc03587; lower panel: 
ugc06138, ugclllSl, ugcll651. 



Finally, we fit point-spread function (PSF) convolved 
King (1966) models to all detected sources using the code 
described in Jordan et al. (2005). The code returns the 
best-fit concentration c and half-light radius r/, for each 
source, but we use only the latter as the former is poorly 
constrained. The half-light radii are useful in selecting 
bona-fide GC candidates from the full list of detections 
(see below). 

3. SURFACE BRIGHTNESS PROFILES AND GALAXY 
CLASSIFICATION 

In order to obtain an improved classification of the tar- 
get galaxies based on the central surface brightness we 
used our observations to derive surface brightness pro- 
files. We then used double Sersic model fits to describe 
the profiles and classify them according to their derived 
values of /Urj,disk) the central surface brightness of the 
large scale (disk) component of the profile. In this section 
we describe these measurements and our new classifica- 
tions. 

3.1. Surface Photometry 

An isophotal analysis of the light profile of the target 
galaxies was performed running the IRAF task "ellipse" 
which is based on the algorithm of Jedrzejewski (1987). 
This task was run on a version of our ACS images where 
bright sources such as foreground stars were masked. We 
used fixed position angles and ellipticities, which were es- 
timated after some experimentation with "ellipse", and 
allowed 20 pixel as the maximum wander between succes- 
sive isophote centers. Using this information we obtained 
an azimuthally averaged intensity profile for each galaxy. 
The "sky" brightness was determined by measuring the 
emission level in five boxes of 200 x 200 pixel 2 located 



in regions of the image that appear to be free of galaxy 
light. The average of the mean counts in the five boxes 
was adopted as the background value for every image and 
the standard deviation between these values was used to 
estimate the errors in the sky determination. 

The measured intensity profiles were converted to AB 
magnitudes per square arcsecond using the photometric 
zeropoints of Sirianni et al. (2005) and a pixel size of 
0'.'05. The profiles were corrected for foreground extinc- 
tion by using E(B — V) values from Schlegel et al. (1998) 
and the extinction ratios listed in Sirianni et al. (2005). 
The resulting surface brightness profiles in the g- and 
i-bands are shown in Figure [21 

3.2. Model Fits and Galaxy Classification 

In order to estimate the surface brightness of the disk 
A*o,disk of our target galaxies we decomposed them by 
fitting two Sersic (1968; see Graham & Driver 2005 for a 
review) profiles, one to describe a bulge component and 
the other to describe a larger scale disk component. A 
Sersic profile is given by 

/(r)=/ e exph6„((r/r e ) 1 /"-l)] (1) 

where n is the Sersic index, I e the intensity of the profile 
at the effective radius r e and b n ensures that r e contains 
half of the integrated light and can be approximated by 
b n w 1.992n - 0.3271 (see Graham & Driver 2005 and 
references therein). 

It is usual to model disks as pure exponential profiles 
(or Sersic profiles with n = 1) but we found that our pro- 
files are usually not well described by exponential pro- 
files on large scales. This is consistent with the poor 
description afforded by exponential fits to describe the 
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Fig. 2. — Surface brightness profiles in the g- (down) and i-bands (up) for the target galaxies. The red-solid line shows the best Sersic 
profile fit to the large (disk) component and the dashed line shows the corresponding inner (bulge) Sersic component, both for the g-band 
profile. The two horizontal lines indicate the limits for normal to intermediate and intermediate to low surface brightness and are located 
at /xb=22 mag/arcsec 2 and p,B=23 mag/arcsec 2 and translated into the g-band by assuming (g — B) = —0.34 (Fukugita et al. 1995). 



large scales of a sample of 21 late- type galaxies with a 
range of surface brightness (Galaz et al. 2006; their on- 
line Figures 2 and 3). We have therefore chosen to use 
the more general Sersic profile to describe the disk com- 
ponent as this can properly describe the data. In all 
cases the best-fit Sersic component describing the disk 
satisfies n < 1. We carried out the fits applying a x 2 
minimization using procedures from CERN's Minuit li- 
brary. The minimization procedure used consists of a 
Simplex minimization algorithm (Nelder & Mead 1965) 
followed by a variable metric method with inexact line 
search (MIGRAD) to refine the minima found by the 
Simplex algorithm. Following Byun et al. (1996) and 
Ferrarese et al. (2006; their equation 14) we weighted all 
points equally in a fractional sense. Our best fit double 
Sersic models are shown in Figure [2] and they reproduce 
the data very closely. The external Sersic component is 
taken to be the disk and to infer the value of /io,disk in the 
g-band. In order to estimate /io,disk in the B-band, we 
assumed a typical color (g — B) = —0.34 mag for an Sbc 
galaxy (Fukugita et al. 1995). The values we inferred for 
the B-band value of no disk are then 21.86 (ugc00477), 
21.43 (ugcO.3459), 22.63 (ugc03587), 23.74 (ugc06138), 
21.62 (ugclllSl) and 21.75 {ugcll651). We used these 
values to re-classify our target galaxies as described be- 
low. 

There is no clear definition of the limiting surface 
brightness separating low from high surface brightness 
galaxies. The convention has been evolving through the 



years as much as the limiting magnitude of the available 
galaxy surveys did. The most recent works on low sur- 
face brightness galaxies are mostly restricted to galaxies 
with central surface brightness of /its > 23 mag/arcsec 2 
(e.g., Zackrisson et al. 2005), but several different defini- 
tions are also used, one of the most common ones being a 
blue central disk surface brightness value of /xo,disk <^ 22 
mag/arcsec 2 (e.g. Boissier et al. 2003). As shown in Fig- 
ure [2] the observed values of /Xo,disk imply that our sam- 
ple covers a large range of central surface brightness. We 
have used the model fits described above to re-classify our 
target galaxies following McGaugh et al. (1995), which 
defines LSBs as those objects showing a central disk sur- 
face brightness in the B band of /x 0i( jisc > 23 mag/arcsec 2 
and "intermediate surface brightness" those galaxies sat- 
isfying 22 mag/arcsec 2 < /io,disc < 23 mag/arcsec 2 . Us- 
ing these definitions we classify ugc06138 as an LSB and 
ugc03587 as an intermediate surface brightness galaxy. 
The other 4 galaxies in our sample are "normal" (or high) 
surface brightness galaxies, despite having low mean sur- 
face brightness. 

The sample of Bothun et al. (1985) from which our 
sample was drawn is selected on mean surface brightness, 
which usually translates into ^io,disc > 23 mag/arcsec 2 
but with objects with values of /xo,disc as high as 22 
mag/arcsec 2 (O'Neil et al. 2004). Our findings are 
roughly consistent with this view when taking into ac- 
count that the four galaxies that fall into the high surface 
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brightness category suffer from a high estimated extinc- 
tion in their field of view. 

4. GLOBULAR CLUSTER CANDIDATE SELECTION 
CRITERIA 

The SExtractor output catalogs of objects in each field 
are strongly contaminated by both foreground Galactic 
stars and background galaxies. Therefore, we need to 
define criteria in order to identify potential globular clus- 
ters from the full list of detections. We selected our GC 
candidates as those objects that fulfill the following con- 
ditions: 

Magnitudes. The GC luminosity function (GCLF; the 
distribution of GCs per unit magnitude) is observed to 
have a nearly gaussian shape with a roughly constant 
peak at a magnitude of My ~ —7.4 mag (Harris 2001), 
which corresponds to M g ~ — 7.1 mag according to stel- 
lar population models for a single-burst population of 10 
Gyr with [Fe/H]= -1.35 (Maraston 2005). For a galaxy 
having a mass like the Milky Way as our targets, the typ- 
ical dispersion a in a Gaussian description of the GCLF 
is a ~ 1.2 mag (Jordan et al. 2006; 2007b). Therefore, 
we adopted an upper luminosity cut at 2.5a above the 
expected GCLF turnover in the (/-band, which translates 
into the condition M g > — 10 mag. 

At low GC luminosities we impose a cut on apparent 
magnitude driven by the completeness limit of the ob- 
servations. The completeness function depend both in 
the magnitude of the objects and their position on the 
galaxy. To account for these effects we added to the 
galaxy fields simulated GCs with a given magnitude and 
a typical size of 3 pc, randomly distributed along the im- 
age. The photometry was then performed in the same 
fashion as for the real observations. After that the num- 
ber of fake added sources recovered was computed. The 
iteration of this process in steps of 0.01 magnitudes at 
least 50 times per field of view allowed us to determine 
the 100% completeness limit of our observation with 10% 
confidence level. These limits correspond to 26.3 mag 
and 25.5 mag in the g- and i-bands respectively and were 
used to fix the low magnitudes cut of our sample selec- 
tion. 

Color. A broad color range, 0.4 < (g — i)q < 2.0, was 
used in order to include GCs with metallicities satisfying 
-2.26 < [Fe/H] < 0.35 for an assumed old r = 10 Gyr 
age and to allow some intrinsic reddening in the target 
galaxies. This color cut corresponds to 0.42 < (B — V) < 
1.43, a range that includes 84% of the galactic GCs with 
available (B — V) color information in the McMaster 
catalog of Milky Way GCs (Harris 1996). 

Morphology. Observations in the Milky Way show that 
GCs are fairly spherical systems (White & Shawl 1987). 
The most extremely flattened Galactic case, NGC 6273, 
has an ellipticity of 0.27, which corresponds to an axial 
ratio of 1.33. Based on these observations we included in 
our catalog only objects whose SExtractor's measured 
axial ratio was a/b < 1.5, leaving an wide range for 
more flattened objects but discarding mostly background 
galaxies. We note that due to the fact that GCs in our 
target are marginally resolved, the adopted condition is 
very conservative as the observed axial ratios are strongly 
affected by the PSF and therefore high axial ratios will 
tend to be circularized. 

We also included a second morphological selection cri- 



teria based on SExtractor's "CLASS_STAR" parameter, 
which on the basis of a trained neural network assigns 
to each source a value between (galaxy) and 1 (star) 
(see Bertin & Arnouts 1996 for details). Objects with 
CLASS-STAR < 0.1 were excluded from the final cata- 
log. 

Size. Using the measured half-light radii we exclude 
from the final catalog all objects that have a half light 
radii in the i band r^i < 0.25 pix. This rejects sources 
that are consistent with being point sources, i.e. mostly 
stars, and therefore eliminates most of the foreground 
contamination in our sample. For the most distant galax- 
ies (D ~ 35 kpc) this cut corresponds to 77^ ~ 2.2 pc 
which would exclude ~ 25% of the total globular clus- 
ters population of the Milky Way if observed at those 
distances because of being more compact than our limit. 
In the case of the closest galaxy (D = 16.83 kpc) this 
value goes down to 1 pc which would exclude just a few 
(~ 5) Galactic GCs. 

We adopted a physical upper limit to the radius of 
Th.i < 10pc, since just ^5% of the Milky Way globular 
cluster systems have r^,i higher than this value. The up- 
per cut in 77, implies that our study does not probe for 
the existence in our sample galaxies of objects such as dif- 
fuse star clusters or faint fuzzies (Brodie & Larsen 2003; 
Peng et al. 2006) and ultra-compact dwarfs or dwarf- 
globular transition objects (Hilker et al. 1999; Drinkwa- 
ter et al. 2000; Hasegan et al. 2005). 

As a final requirement on the structural parameters, we 
eliminated objects whose fractional uncertainties Sr^ = 
<Jr h / r hi where a rh is the uncertainty in the measured ru, 
satisfy Srh > 0.3 in any of the bands. A large error in 
the measured is an indication that the observed light 
distribution is not properly fit by a King (1966) model, 
and as a consequence it is likely that the source is a 
background galaxy instead. 

Figure [3] shows the location in a color magnitude dia- 
gram of objects classified as globular cluster candidates 
according to the criteria described above (red crosses). 
The small background dots in the same diagram repre- 
sent all additional sources detected in the field of view of 
each galaxy. 

4.1. Control Fields 

Even though we expect to obtain a fairly clean cat- 
alog of GC candidates after applying the selection cri- 
teria described above, some contamination will remain 
in the sample. In order to estimate the expected level 
of remaining contamination in the object sample of each 
galaxy, we searched the HST/ACS archive for blank fields 
at least as deep as our observations, and whose galactic 
coordinates would make them good control samples, ei- 
ther because they are located close to one of our fields, 
or in a symmetric position with respect to the Galac- 
tic center. These conditions were chosen such that the 
properties of the galactic foreground in the control fields 
were similar to our observed frames, but note that given 
that we can eliminate most stars using our half-light radii 
measurements a very accurate match is not crucial. Two 
of the galaxies, ugcOOJ^ll and ugc06138, are located at 
very high galactic latitude, therefore any high galactic 
latitude blank field is suitable to reproduce the contami- 
nation as in this case the amount of foreground Galactic 
stars is very low. For these galaxies we were able to ob- 
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Fig. 3. — Color-magnitude diagrams of all the sources detected in the field of view of each galaxy (small background dots). The red 
crosses correspond to those sources classified as globular cluster candidates. 



TABLE 3 
Control Fields 



ID 


I 


b 


Prog. ID 


E(B-V) a 


C01 


145.58 


-38.65 


9488 


0.155 


C02 


115.01 


46.68 


9488 


0.018 


C03 


99.41 


-47.15 


9488 


0.033 


C04 


289.08 


62.25 


9488 


0.024 


C05 


251.33 


-41.44 


9488 


0.008 


C06 


92.66 


46.37 


9488 


0.016 


C07 


279.93 


-19.99 


9575 


0.130 


C08 


179.06 


-19.93 


9488 


0.630 


C09 


105.10 


7.07 


9488 


0.824 


CfO 


24.48 


12.67 


9488 


0.478 



a Foreground Galactic reddening in the direction of each field of view 
from Schlegcl et al. 1998. 

tain six control fields which are denoted C01 to C06 in 
Table [3J The contamination for these two galaxies was 
estimated as the average expected number of contami- 
nating objects in the six control fields. 

The other four galaxies in the sample have low galactic 
latitude, making the selection of control fields harder as 
not many suitable observations have been carried out in 
such regions with the filter set used in this program. In 
these cases we had to restrict ourselves to finding a con- 
trol field as close as possible to our target, which meant 
using just one or two control fields per galaxy. Column 2 
in Table 0] indicates the ones we consider to be the best 
control field for a given galaxy according to the labels 



used in Table [3J 

All the chosen control fields were observed as part of 
two ACS Pure Parallel programs, GO-9488 and GO- 
9575, and they have exposure times 1600 sec < t < 
1800 sec in both bands. These images were reduced and 
the photometry was performed using exactly the same 
procedures followed for our targets. In order to account 
for the influence of the light of each target galaxy on the 
detections over our frames we follow a procedure similar 
to that described in §2.2 of Peng et al. (2006a). Namely, 
for each galaxy we created "custom" control samples by 
re-doing the detections as if a given galaxy was in front 
of our blank control fields. To do this, we used the noise 
images created for each galaxy in the data reduction pro- 
cess when running the detection procedure on the con- 
trol fields. In this way we can reproduce for each control 
field the catalog of objects that we would have detected 
if there had been a particular galaxy of our sample in the 
field. 

The same selection criteria described above were ap- 
plied to the final SExtractor catalogs of the control fields 
in order to obtain a final catalog of the expected con- 
taminating objects. As we aimed at reproducing the ex- 
pected contamination in the case the control field was 
in the same position as our target fields, the apparent 
magnitude for the control objects was computed as: 

m = m obs - A cf + Agf (2) 
where m \, s is the observed magnitude in a given band 



Villegas et al. 



and A c f and A g f are the extinction corrections for the 
control field and the galaxy field in the same band re- 
spectively. 

5. TOTAL NUMBER OF GLOBULAR CLUSTERS 

Once the selection criteria have been applied to all the 
objects in each field we have the total number of ob- 
served GC candidates. In order to have an accurate 
estimate of the total population of GCs in our target 
galaxies, we need to consider both the expected number 
of contaminants and the fraction of objects we might be 
losing due to observational and selection effects. Below 
we describe our adopted procedures to account for these 
observational effects. 

5.1. Contamination 

There are mainly two different kinds of contamination 
in our sample. One is external in nature and is composed 
mostly of background galaxies and any foreground star 
whose inaccurate Th measurement might have allowed it 
to pass our rejection of unresolved sources. This exter- 
nal contamination can be corrected by using the control 
fields customized to the field of each target galaxy. An- 
other set of contaminants is of internal nature, and is 
composed of objects in each target galaxy which given 
our adopted selection criteria in morphology, color, mag- 
nitude and size are almost certainly reddened star form- 
ing complexes. 

The typical size for an ultra compact H II region is 
r ^ 1-6 pc (Churchwell 2002). Objects with these sizes 
might be present in our sample in the three closest galax- 
ies, and this could constitute a source of contamination, 
as some objects with those sizes are indeed detected (see 
Figure [5] and discussion below). The other main ex- 
pected contribution from non-GC sources in the target 
galaxies to our candidate sample will be due to reddened 
young star clusters lying in the spiral arms of our tar- 
gets. Young star clusters typically show sizes between 
1-20 pc, with typical radii of r ~ 3 pc (e.g., Scheep- 
maker et al. 2007), magnitudes satisfying My > —11.5 
and colors in the range —0.2 < (B — V) < 0.4 (Larsen 
1999), roughly corresponding to —0.57 < (g — i) < 0.36. 
Therefore, even a small amount of reddening may allow 
young star clusters to be included in our sample of (old) 
GC candidates. Visual inspection of our images shows 
that most of our GC candidates are not located in obvi- 
ous star forming regions, even though a small amount of 
contamination is still possible. 

We show in SJ5] that the properties (luminosity func- 
tion, color and size distribution) of the combined sample 
of GC candidates in our sample of LSB galaxies are con- 
sistent with the expected properties of old GCs in HSB 
galaxies (see below). Therefore, we are confident that 
our sample of GC candidates is not dominated by young 
star clusters. Of course, while the total sample is likely 
dominated by old GCs, the nature of any given object 
can only be assessed by spectroscopic follow-up. 

5.2. Incompleteness 

Due to observational effects and the various selection 
criteria we are using, we do not expect to be able to 
observe the complete GC system of each galaxy. There- 
fore, we need to estimate the fraction of the GC system 
we might be missing and use this information to estimate 



the total size of the GC systems. As the Milky Way has 
the best studied GC system and its size and morphol- 
ogy are comparable to those of our target galaxies, we 
used the observed properties of its GC system in order 
to estimate the total number of clusters in each galaxy 
through a direct comparison. In order to do so, we follow 
a procedure similar to the one described by Kissler-Patig 
et al. (1999). Based on the areal coverage of the WFC 
we create a spatial "mask" for each galaxy and we apply 
this mask and our additional selection criteria to a pro- 
jection of the Galactic GC system as described below. 
The properties of the Galactic GC systems are obtained 
form the McMaster catalog (Harris 1996). 

We apply our selection criteria to the two-dimensional 
spatial distribution obtained projecting the Galactocen- 
tric coordinates (X, Y, Z) of the Galactic GC system into 
the Y-Z plane and rejecting sources that lie outside the 
spatial mask of each galaxy and that do not satisfy the 
selection criteria presented in Sj4j The total number of 
GCs for a given galaxy, Nqc, is then given by 

N GC = N M w^ £ ^, (3) 

where Nqc is the number of GC candidates in each tar- 
get galaxy left after subtracting from the observed num- 
ber of GC candidates the external contamination, iV mas k 
is the number of Galactic clusters detected inside the 
mask and Nmw is t ne total number of clusters in the 
Milky Way which we assume to be 180±20 (Ashman & 
Zepf 1998). In order to apply to the Milky Way GCs the 
same selection criteria we apply to our target galaxies, 
the cuts in magnitude were applied to absolute magni- 
tude of the Galactic clusters that were reddened as ap- 
propriate for each of our target galaxies, and (B — V) 
colors were transformed to (g — i) using a linear conver- 
sion derived from the Maraston (2004) models 8 . 

This direct comparison with the GC system of our 
Galaxy relies on several assumptions which we now dis- 
cuss. 

5.2.1. The Galactic GC sample 

The McMaster Catalog lists a total of 150 globular 
clusters, all of them with Galactocentric coordinates. 
There are four GCs without photometric information, so 
they are excluded from the comparison sample and we 
start working with 146 object, accounting the 4 missing 
clusters as part of the uncertainty in the total number 
of GCs in the Galaxy. From these 146 GCs, 29 have no 
color information and, from the remaining sample, 17 do 
not have measured ellipticities. Therefore, the compari- 
son sample is limited to 100 Galactic GCs. 

The exclusion of these 46 clusters does not represent 
a significant bias in the case of ugc00477, ugc03459, 
ugc06138 and ugclllSl, because at the distances and 
reddening factors of these galaxies those objects would 
be rejected anyway either by the faint magnitude cut or 
by the color/ellipticity selection, i.e. they are rejected 
without the need to use the missing information. For 
ugc03587 and ugcll651 we face a different situation, be- 
cause after applying all the selection criteria, except the 
one that uses the missing quantity to these 46 objects, 14 

8 The transformation used was [g — i) = 1.57(-B — V) + 0.26. 
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and 5 objects are still left for those two galaxies respec- 
tively. These numbers constitute an additional source of 
uncertainty for these galaxies. 

Following Ashman & Zepf (1998), we assume that the 
Milky Way has a total of 180±20 GCs, where the un- 
observed clusters are assumed to be hidden behind the 
Galactic bulge. We assume in what follows that in our 
target galaxies the fraction of obscured GCs is the same 
on all galaxies and is given by the Galactic estimate of 
this fraction. 

5.2.2. The Magnitude Distribution of the Milky Way GCs 

When setting up the selection criteria and when com- 
paring the GC population of the target galaxies with the 
Milky Way's GC system we are implicitly assuming that 
the luminosity function of both systems is similar. It has 
been recently shown that the GCLF depends systemati- 
cally on the galaxy mass (Jordan et al. 2006; 2007b), but 
given that all of our target galaxies have similar masses 
between them and with the MW we can safely assume 
that the GCLFs of our target galaxies are roughly similar 
than that of the MW. In order to make a consistent com- 
parison between the systems, we need to make sure to 
consider the same fraction of the GCLF in the Milky Way 
and any sample galaxy. In order to do that the cut at 
faint magnitudes was applied over the (un-de-reddened) 
absolute magnitudes of the Galactic clusters and then 
dimmed by the distance and absorption coefficients cor- 
responding to each galaxy. At bright magnitudes the cut 
was applied at My ~ —10.3. 

5.2.3. The Spatial Distribution of the Milky Way GCs 

Another implicit assumption done when using the spa- 
tial mask is that the spatial distribution of the Milky 
Way GC system is similar to that of the GC systems 
of our target galaxies. Little is known about the spa- 
tial distribution of GCs in spiral galaxies. In general, 
the spatial distributions seem to be properly fitted by a 
projected power law with indices ranging from -2 to -2.5, 
but some differences exist between Milky Way-like galax- 
ies with a more concentrated distribution, and M31-like 
galaxies where the GC system is more extended than the 
halo light distribution (see Ashman & Zepf 1998 for a 
review). Harris (1986) and Kaisler et al. (1996) have 
shown that more luminous galaxies have shallower radial 
distributions. 

If the sample galaxies had much steeper profiles as 
compared to the one observed in the Galaxy we would be 
overestimating the number of clusters that are left out of 
the mask, but this would not greatly affect the final mea- 
surement, because even for the galaxy with the smaller 
areal coverage (ugc0358T) the mask includes more than 
half of the the Galactic globular clusters, when assuming 
it has a similar profile. If the spatial distribution of GCs 
in the galaxies happen to be shallower than the observed 
in the Milky Way we will be accordingly underestimating 
the total number of clusters of the system. 

5.2.4. Distances 

The uncertainty in the distances to the galaxies is likely 
to be the most important systematic error in this study, 
because distances are involved in the absolute magnitude 
and size selection and they are also crucial to accurately 



reproduce the field of view for each galaxy when creat- 
ing the comparison mask. The distances we used are 
corrected for the Virgo infall, but still have the intrinsic 
error due to the peculiar velocity of field galaxies and the 
uncertainty in the value of Hq, which we assume to be 73 
km/s/Mpc. Estimates of the peculiar velocity for field 
galaxies depend on the method and the sample used, but 
they are thought not to exceed ~ 100 km sec -1 (see, e.g., 
Davis et al. 1997). 

As no better distance estimate of our target galaxies is 
available, we evaluated how our results change under a 
±15% variation in the assumed distances. This exercise 
is then used to account for the uncertainty due to peculiar 
velocities and the error in Hq. 

5.3. The Total Number of GCs 

We determined the total number of clusters for each 
galaxy following the method described above. Table [4] 
shows our results and their uncertainties using our as- 
sumed distances to the galaxies as well as the results 
with these distances varied by ±15%. The uncertain- 
ties reported in Table |4] include the random errors de- 
scribed above estimated assuming Poisson statistics for 
the number of GC candidates, the expected contamina- 
tion in each field and the number of objects inside the 
Milky Way's mask. 

Additionally, several sources of systematic uncertain- 
ties need to be accounted for. To the uncertainty in the 
distances (which introduces significant differences in the 
final number of clusters), we need to add the 14 and 5 
Galactic GCs lacking information in the case of ugc03587 
and ugcll651, respectively. The difference in the final 
GC numbers due to these missing objects should not ex- 
ceed -22/+32 and -6/±6 objects for each of these galax- 
ies, respectively. Furthermore, we assumed a fraction of 
~ 20% of obscured clusters in both galaxies even though 
this may not be the case. Finally, but certainly not less 
important, four galaxies in the sample are located at very 
low galactic latitudes, where we have to work with large 
and uncertain extinction values. Especially critical is the 
case of ugc03459, where the high galactic extinction com- 
bined with the distance to the galaxy led in practice to 
a result consistent with a null detection of cluster candi- 
dates. 

6. GENERAL PROPERTIES OF THE GLOBULAR CLUSTER 

SYSTEMS 

In what follows we present general properties of the 
GC systems detected in our program galaxies. When- 
ever appropriate, we divide the sample according to our 
derived values of /io.disk in order to probe for any effect 
that surface brightness might have on the properties un- 
der study. In all the following sections except in ^6 . 11 
we refer to the properties of all the 206 objects classi- 
fied as globular cluster candidates, without performing a 
statistical subtraction of the expected contamination. 

6.1. Specific Frequencies 

We present in Table [4] the derived specific frequen- 
cies for our target galaxies, which are defined as Sn — 
N GC 10 OA ( M v+i5) (Harris & van den Bergh 1981), where 
My is the total visual magnitude of the galaxy which was 
transformed from the known Mb magnitude assuming a 
typical color for an LSB galaxy of (B — V) = 0.6 mag (see 
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Table 2 in de Blok et al. 1995). The statistical error was 
computed in this case assuming a 0.2 mag uncertainty in 
the magnitude of the galaxy. 

Within the uncertainties the results are consistent with 
the typical Sn for late-type galaxies, i.e. Sn ~ 1- The 
average Sn for the 5 sample galaxies with GC detections 
is Sn = 2.34 ± 0.99, which is slightly higher than the 
average Sn for late-type galaxies. McLaughlin (1999) 
suggested from studies in early-type galaxies that there 
is a universal GC formation efficiency by mass e c i (i.e. the 
fraction of the total baryonic mass that turns into GCs; 
see Blakeslee 1999 for a similar proposal), and inferred a 
value of e c i w 0.25%. The fact that the specific frequen- 
cies we measured are a factor of ~ 2 higher than those 
typical in normal brightness late-type galaxies is consis- 
tent with the idea of a universal GC formation efficiency 
given the fact that galaxies classically classified as hav- 
ing low surface brightness have mass-to-light ratios that 
are « 2 times higher than of higher surface brightness 
galaxies (Zwaan et al. 1995; Impey & Bothun 1997). In 
particular, the specific frequency for our two lowest sur- 
face brightness galaxies (ugc03587 and ugc06138) is ~ 2, 
albeit with large uncertainties. 

6.2. Color Distribution 

The left panel of Figure [4] shows the color distribu- 
tion of the 206 globular cluster candidates detected in 
six galaxies as well as the joint color distribution of 
the two intermediate to low surface brightness galaxies 
(ugc03587 and ugc06138) and the individual distribution 
of ugc06138. The colors in this panel are not corrected 
for internal galaxy reddening. The right panel shows 
the (un-de-reddened) color distribution of the Milky Way 
GCs that were detected inside the comparison masks of 
the six galaxies. In other words, the right panel shows 
how the left panel would look like if the GCs in all the 
galaxies would have color distributions such as the one 
of the Milky Way. 

Comparing the total sample color distribution to the 
one corresponding to the Milky Way, we note that there 
is a higher number of GCs in the bluest bin of our sam- 
ple, probably due to our relaxed color selection crite- 
ria, which can lead to the inclusion of some young star 
clusters in our sample. The results of a Kolmogorov- 
Smirnov test performed in order to compare both dis- 
tributions indicates that the null hypothesis that both 
samples arise from the same parent distribution cannot 
be rejected (p- value of 0.4) when restricting the color to 
satisfy (g — i) > 0.68 mag, a regime where both distribu- 
tions overlap. We conclude then that the GC candidates 
observed in the LSB galaxies have a color distribution 
similar to the one observed in the Milky Way in the re- 
gion where they overlap. An additional group of objects 
sharing the general features of GCs but with bluer col- 
ors is present in the sample galaxies. Whether this color 
is due to low metallicities or young age is something we 
cannot answer with our data but can be addressed with 
follow-up spectroscopy. 

With respect to the color distribution of ugc03587 and 
ugc06138, the data suggest that they roughly follow the 
same distribution as the rest of the galaxies in the sam- 
ple. While no evident differences are apparent as a func- 
tion of central surface brightness, we note that we are 
dealing with small number statistics. 
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Fig. 4. — Left: The dotted line show the color distribution of 
globular cluster candidates in the 6 galaxies, compared to the one 
observed for the 2 intermediate to low surface brightness galaxies 
(ugc03587 and ugc06138) in the solid line, and to the one observed 
for ugc06138 (dashed line) . Right: Color distribution of the Galac- 
tic globular clusters detected under the same observing conditions 
and selection criteria we applied for the six galaxies. 



It has been known for some time that the color dis- 
tributions of early type galaxies are bimodal (see, e.g. , 
Peng et al. 2006b and references therein) . This property 
has also been observed in the GC populations of normal 
spiral galaxies, including the Milky Way and M31. For 
the latter galaxies, the color bimodality has been shown 
to correspond to metallicity bimodality (see e.g., Ashman 
& Zepf 1998, Barmby et al. 2000). The number of GC 
candidates in our galaxies is not high enough to differen- 
tiate between uni- or bimodality, but we note that some 
asymmetry can be observed in the color distributions. 

6.3. Size Distribution 

The half light radii distribution of GC candidates is dis- 
played in Figure [5] for all the galaxies in the sample and 
for the same sub-samples in surface brightness described 
in the previous section. We have overplotted in this figure 
the analytical distribution that Jordan et al. (2005) have 
shown to be a good fit for the size distribution of GCs in 
the Milky Way and in early-type galaxies, convolved with 
the error distribution observed in our sample. We are us- 
ing their eq. 22 with parameters: fi = 0.37, f3\ — 0.117, 
/3 2 = 0.078 and / = 0.7. 

Our sample shows a size distribution which is qualita- 
tively similar to the one observed in early-type galaxies 
and the Milky Way (Jordan et al. 2005). In particular, 
this distribution shows a peak at a half-light radius of 
~ 3 pc, with an extended tail to larger distances. We 
note that the fact that the typical sizes of our GC candi- 
dates are ^ 3 pc is an independent confirmation that our 
assumed distances are reasonable. As shown in Jordan 
et al. (2005) following previous suggestions (e.g., Kundu 
& Whitmore 2001), the typical sizes of GCs can be used 
as a standard ruler and therefore biased distances would 
have returned typical GC sizes in disagreement with the 
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TABLE 4 

Number of globular clusters and upper limit for the specific frequency per galaxy." 



15% closer 15% farther 



Galaxy 


Control Field(s) 


Ncc 


Nmw 


N G c 


Sn 


N G c 


Sn 


N G c 




(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(7) 


(8) 


(9) 


(10) 


UGC 00477 


C01-C06 


42 (49) 


31 


244 ± 80 


3.10 ± 1.17 


172 ± 54 


2.27 ± 0.83 


263 ± 94 


3.46 ± 1.39 


UGC 03459 


C08,C09 


(2) 


6 


- ± - 


- ± — 


- ± - 


- ± — 


- ± - 


- ± - 


UGC 03587 


C08,C09 


57 (61) 


57 


180 ± 43 


2.37 ± 0.72 


177 ± 43 


2.33 ± 0.71 


162 ± 38 


2.13 ± 0.63 


UGC 06138 


C01-C06 


13 (20) 


31 


75 ± 51 


2.02 ± 1.41 


46 ± 36 


0.86 ± 0.70 


94 ± 58 


1.75 ± 1.12 


UGC 11131 


C10 


19 (25) 


19 


180 ± 81 


1.53 ± 0.75 


173 ± 67 


1.46 ± 0.63 


219 ±112 


1.85 ± 1.01 


UGC 11651 


C07,C09 


42 (49) 


38 


199 ± 57 


2.72 ± 0.93 


143 ± 41 


1.96 ± 0.66 


201 ± 59 


2.74 ± 0.95 



a (1) Galaxy ID. (2) Selected control fields as named on Table 2. (3) Number of detected cluster candidates after subtraction of the expected 
contamination. The total number of GC candidates detected is indicated in brackets. (4) Number of MW globular clusters that would be detected 
under the same observational conditions. (5) Expected total number of clusters after comparison with the MW system. (6) Specific frequency 
derived from column 6. (7-8) and (9-10) The same as (5-6) but for a 15% smaller and 15% larger galaxy distance respectively. 
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Fig. 5. — Size distribution of globular cluster candidates de- 
tected in the 6 galaxies (dotted line) compared to the distribution 
observed for the 2 intermediate to low surface brightness galaxies 
(ugc03587 and ugc06138) in the solid line, and to the one observed 
for ugc06138 (dashed line). The dot-dashed line shows a fit for the 
size distribution of GCs in the Milky Way convolved with the error 
distribution observed in our sample (see text). 



expected values. 

In the small half light radius regime we observe a group 
of very compact objects (77, <1.5) that do not seem to fit 
properly the Galactic size distribution. This population 
overlaps with the one responsible in creating an excess 
at the blue region of the color distribution, which in- 
dicates that they are most likely internal contamination 
due to young compact clusters in the disk of the galaxies, 
whose color and size distributions happen to overlap the 
ones for globular clusters in these regions (Larsen 1999, 
Scheepmaker et al. 2007). 

6.4. Luminosity Function 

Figure [6] shows the luminosity distribution of the glob- 
ular cluster candidates in all 6 galaxies in both bands. 
This figure also shows the luminosity distribution of the 
sub-sample of intermediate to low luminosity galaxies 
and the individual distribution of ugc06138. The Gaus- 
sian distribution over-plotted in both panels of Figure [5] 
shows the expected globular cluster luminosity function 



for a galaxy of Mb ~ -19, whose parameters correspond 
to (M g , a g )=(-7.2, 1.2) and (Mi, cr i )=(-8.1, 1.2), accord- 
ing to Jordan et al. (2007b). These distributions were 
normalized to twice the number of cluster candidates 
brighter than the peak (turnover) of the luminosity func- 
tion. By normalizing in this fashion we circumvent the 
effects of incompleteness present at magnitudes fainter 
than the turnover. 

The luminosity distribution of all the cluster candi- 
dates gives the appearance of being brighter than ex- 
pected in both bands. This is due to different fore- 
ground absorption conditions between the sample galax- 
ies, which do not allow to reach equally deep magnitudes 
for all of them. In fact, ugc03587 is the only case where 
the sample extends ^2.5<t down the expected peak of 
the distribution. The samples of ugc00477, ugc06138 
and ugc 11651 reach at least the turnover, and in the 
cases of ugc03459 and ugclllSl we are sampling just the 
most luminous objects. Considering this caveat, the ob- 
served distributions are consistent with the expectations 
from the luminosity function of the MW and early-type 
galaxies (e.g. Jordan et al. 2007b). 

It is also worth noticing that the luminosity distribu- 
tions of clusters hosted by ugc03587 and ugc06138 are 
consistent with those of the GCs hosted by high surface 
brightness galaxies in our sample. This implies, under 
the standard picture of the dynamical formation of the 
GC luminosity function (e.g., Fall & Zhang 2001), that 
two-body relaxation and evaporation have acted with 
similar effectiveness in both kind of systems. Given that 
the evaporation rate can be argued to depend mainly on 
the GC half-mass density (e.g., McLaughlin & Fall 2007) 
we expect the extent of GC dynamical evolution to be 
similar across our sample given that the luminosity and 
size distributions show no systematic variations, imply- 
ing the same for the GC half-mass densities. 

7. SUMMARY AND CONCLUSIONS 

We have studied the GC systems of a sample of six 
massive late-type galaxies with a range of surface bright- 
ness, including two galaxies with blue central disk surface 
brightness /Jo, disc > 22 mag/arcsec 2 . 

Through direct comparison with the Milky Way, we 
estimated the total number of globular clusters and us- 
ing these values we have derived specific frequencies for 
each galaxy. The results obtained do not show significant 
differences between galaxies in different ranges of surface 
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Fig. 6. — Luminosity distribution of the globular cluster candi- 
dates following the same line-style used in the previous 2 figures. 
The dot-dashed line shows the expected GCLF for a galaxy with 
Mg ~ 19, a Gaussian distribution with parameters (M g , <r s )=(- 
7.2, 1.2) and (A?,, cr i )=(-8.1, 1.2) (Jordan et al. 2007b). 

brightness and also indicate that all of them have formed 
globular clusters with similar efficiency in the range ex- 
pected for late-type galaxies. 

A general analysis of the properties of the 206 clusters 
candidates detected in the 6 sample galaxies has shown 
that they all show similar characteristics as those ob- 
served in previously studied GC systems. Their lumi- 
nosity function turns over close to the expected value of 
M g ~ —7.2 and their size distribution peaks at ~3 pc, 
just like in the case of the Milky Way and many other 
spiral and elliptical galaxies. Their color distribution is 
similar to the one expected if the Galactic globular clus- 
ters were observed under the same conditions as the stud- 
ied galaxies. All these properties do not show significant 
evolution with the surface brightness of the host galaxy. 

When the properties of two intermediate-to-low sur- 
face brightness galaxies are compared to the high sur- 
face brightness galaxies in the sample they do not show 
particular features. In all the properties we probed the 
globular cluster candidates in ugc03587 and ugc06138 
show no significant differences with respect to the other 
galaxies even though they have a much lower central 
surface brightness. This implies that the physical con- 
ditions present during the formation of GCs in massive 
LSB galaxies were not significantly different to those of 
HSB galaxies. 

Our results show that the GC formation process in LSB 
galaxies proceeded in a way that, at least with the avail- 
able information, is indistinguishable from the processes 
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observed in normal (i.e high surface) brightness galaxies. 
Our data suggests a close similarity in the early assembly 
processes of both types of galaxies, which implies that de- 
spite a very quiescent star formation history thereafter, 
the initial period of assembly must have been intense 
enough to allow the formation of massive star clusters. 
Improved constraints can be obtained after studying the 
spectra of the globular clusters in order to date their for- 
mation time and measure the metal abundances of the 
environment in which they where formed. 

We have found the presence of some "contamination" 
in the sample of cluster candidates, coming from a group 
of blue and small (r^ <1.5) objects. We are unable to 
determine whether the blue color is due to low metal- 
licity or young ages. The first possibility would mean 
that these galaxies were able to form lower metallicity 
objects than the ones observed in the Milky Way, while 
the second one would imply that ongoing clustered star 
formation in the observed galaxies (and particularly in 
ugc03587 which is the main contributor to this popula- 
tion) is driven in a fashion similar to that observed in 
normal brightness galaxies, despite the low-star forma- 
tion rates expected based on the Kennicutt law (Kenni- 
cutt 1989). 

It is worth mentioning that all the 6 galaxies included 
in our sample were originally cataloged as "low surface 
brightness galaxies" , but after a detailed decomposition 
of their surface brightness profiles we have found just two 
of them to be consistent with this classification. This 
stresses the importance of performing a surface bright- 
ness decomposition and proper absorption corrections 
when classifying objects as LSB galaxies. 

The results we have presented here are broadly con- 
sistent with the scenario proposed by van den Hoek et 
al. (2000) in which LSB galaxies roughly follow the same 
evolutionary history as HSB galaxies except at a much 
lower rate. We have show that the early formation pro- 
cess must have been intense enough in both types of 
galaxies to allow the formation of massive star clusters, 
even though the ensuing evolution proceeds at very dif- 
ferent rates. 



Support for program GO-10550 was provided through 
a grant from the Space Telescope Science Institute, which 
is operated by the Association of Universities for Re- 
search in Astronomy, Inc., under NASA contract NAS5- 
26555. This research has made use of the NASA/IPAC 
Extragalactic Database (NED) which is operated by the 
Jet Propulsion Laboratory, California Institute of Tech- 
nology, under contract with the National Aeronautics 
and Space Administration. 

Facility: HST (ACS/WFC) 



Boissier, S., Monnier Ragaigne, D., Prantzos, N., van Driel, W., 
Balkowski, C. & O'Neil, K., 2003, MNRAS, 343, 653 

Bothun, G.D., Beers, T.C., Mould, J.R. & Huchra, J. P., 1985, AJ, 
90, 2487 

Bothun, G.D., Schombert, J.M., Impey, CD., Sprayberry, D. & 
McGaugh, S.S., 1993, AJ, 106, 530 



LSB Galaxies and their GCs 



13 



Byun, Y.I., Grillmair, C.J., Faber, S.M., Ajhar, E.A., Dressier, A., 
Kormendy, J., Lauer, T.R., Richstone, D.& Tremaine, S., 1996, 
AJ, 111, 1889 

Churchwell, E., 2002, ARA&A40, 27 

Cote, P., Blakeslee, J. P., Ferrarese, L., Jordan, A., Mei, S., Merritt, 
D., Milosavljevic, M., Peng, E.W., Tonry, J.L., & West, M.J., 
2004, ApJS, 153, 223 

Dalcanton, J. J., Spergel, D.N. & Summers, F.J., 1997, ApJ, 482, 
659 

Davis, M., Miller, A & White S.D.M., 1997,ApJ, 490, 63 
Disney, M., 1976, Nature, 263, 573 

de Blok, W.J.G., van der Hulst, J.M. & Bothun, G.D., 1995, 

MNRAS, 274, 235 
de Blok, W.J.G., van der Hulst, J.M., 1998, A&A, 335, 421 
Ferrarese, L., Cote, P., Jordan, A., Peng, E.W., Blakeslee, J. P., 

Piatek, S., Mei, S., Merritt, D., Milosavljevic, M., Tonry, J.L., & 

West, M.J. 2006, ApJS, 164, 334 
Freeman, K.C., 1970, ApJ, 160,811 

Fukugita, M., Shimasaku, K. & Ichikawa, T. 1995, PASP, 107, 945 
Galaz, G., Villalobos, A., Infante, L. & Donzelli, C, 2006, AJ, 131, 
2035 

Graham, AW. & Driver, S.P., 2005, Publ. Astron. Soc. Australia, 
22, 118 

Harris, W. E. & van den Bergh, S., 1981, AJ, 429, 177 
Harris, W.E., 1991, ARA&A, 29, 543 
Harris, W.E., 1996, AJ, 112, 1487 

Harris, W.E., 2001, in Star Clusters, ed. L. Labhardt & B. Binggcli 

(Berlin: Springer), 223 
Impey, C, & Bothun, G. 1997, ARA&A, 35, 267 
Jedrzejewski, R. I., 1987, MNRAS, 226, 747 

Jimenez, R., Padoan, P., Matteucci, F., & Heavens, A. 1998, 

MNRAS, 299, 123 
Jordan, A., Blakeslee, J. P., Peng, E.W., Mei, S., Cote, P., Ferrarese, 

L., Tonry, J.L., Merritt, D., Milosavljevic, M., & West, M.J. 2004, 

ApJS, 154, 509 

Jordan, A., Cote, P., Blakeslee, J. P., Ferrarese, L., McLaughlin, D., 
Mei, S., Peng, E.W., Tonry, J.L., Merritt, D., Milosavljevic, M., 
Sarazin, C.L., Sivakoff, & West, M. 2005, ApJ, 634, 1002 

Jordan, A., McLaughlin, D.E., Cote, P., Ferrarese, L., Peng, E.W., 
Blakeslee, J. P., Mei. S., Villegas, D., Merritt, D., Tonry, J.L., & 
West, M.J. 2006, ApJ, 651, L25. 

Jordan, A., Blakeslee, J. P., Cote, P., Ferrarese, L., Infante, L., Mei, 
S., Merritt, D., Peng, E.W., Tonry, J.L. & West, M.J., 2007a, 
ApJS, 169, 213 

Jordan, A., McLaughlin, D.E., Cote, P., Ferrarese, L., Peng, E.W., 
Mei. S., Villegas, D., Merritt, D., Tonry, J.L., & West, M.J. 
2007b, ApJS, in press. 

Kaisler, D., Harris, W.E., Crabtree, D.R. & Richer, H.B., 1996, 
AJ, 111, 2224 

Kennicutt, R.C., 1989, ApJ, 344, 685 

King, I.R., 1966, AJ, 71, 64 

Kissler-Patig, M., Ashman, K.M., Zepf, S.E. & Freeman, K.C., 
1999, AJ, 118, 197 



Koekemoer, A.M., Fruchter, A.S., Hook, R.N. & Hack, W., 
2002, in The 2002 HST Calibration Workshop. Ed. S. Arribas, 

A. Koekemoer & B. Whitmore, (Baltimore, STScI), 337 
undu, A. & Whitmore, B.C., 2001, 121, 2950 

Larsen, S.S. 1999, A&AS, 139, 393 
McLaughlin, D.E., 1999, AJ, 117, 2398 
Maraston, C, 1998, MNRAS, 300, 872 
Maraston, C, 2005, MNRAS, 362, 799 

McGaugh, S.S., Schombert, J.M. & Bothun, G.D., 1995, AJ, 109, 
2019 

McGaugh, S.S., Bothun, CD. & Schombert, J.M., 1995, AJ, 110, 
573 

Mo, H.J., McGaugh, S.S. & Bothun, G.D., 1994, MNRAS, 267, 129 
Mould, J.R., Huchra, J. P., Freedman, W.L., Kennicutt, R.C., 
Ferrarese, L., Ford, H.C., Gibson, B.K., Graham, J. A., Hughes, 
S.M.G., Illingworth, G.D., Kelson, D.D., Macri, L.M., Madore, 

B. F., Sakai, S., Sebo, K.M., Silbermann, N.A. & Stetson, P.B., 
2000, ApJ, 529, 786 

Nilson, P., 1973, Uppsala General Catalogue of Galaxies, Uppsala 

Astr. Obs. Ann., v. 6 
Nelder, J. A., & Mead, R. 1965, Comput. J., 7, 308 
O'Neil, K., Bothun, G., van Driel, W. & Monnier Ragaigne, D., 

2004, A&A, 428, 823 
Peng, E.W., Jordan, A., Cote, P., Blakeslee, J. P., Ferrarese, L., 

Mei, S., West, M.J., Merritt, D., Milosavljevic, M., & Tonry, 

J.L., 2006a, ApJ, 639, 95 
Peng, E.W., Cote, P., Jordan, A., Blakeslee, J. P., Ferrarese, L., 

Mei, S., West, M.J., Merritt, D., Milosavljevic, M. & Tonry, J.T., 

2006b, ApJ, 639, 838 
Rosenbaum, S.D. & Bomans, D.J., 2004, A&A, 422 
Scheepmaker, R.A., Haas, M.R. , Gieles, M., Bastian, N., 

Larsen, S.S. fc Lamers, H.J.G.L.M., 2007, A&Ain press 

(arXiv: 0704. 3604) 
Schombert, J.M., Bothun, G.D., Impey, CD. & Mundy L.C, 1990, 

AJ, 100, 1523 

Schombert, J.M., Bothun, CD., Schneider, S.E. & McGaugh, S.S., 

1992, AJ, 103, 1107 
Schlegel, D.J., Finkbeiner, D.P. & Davis, M., 1998, ApJ, 500, 525 
Sharina, M.E., Puzia, T.H. & Makarov, D.I., 2005, A&A, 442, 85 
Sirianni, M., et al. 2005, PASP, 117, 1049S 

van den Hoek, L.B., de Blok, W.J.C, van der Hulst, J.M. & de 

Jong, T., 2000, A&A, 357, 397 
van der Hulst, J.M., Skillman, E.D., Smith, T.R., Bothun, CD., 

McGaugh, S.S. & de Blok, W.J.C, 1993, AJ, 106, 548 
West, M.J., Cote, P., Marzke, R.O., & Jordan, A. 2004, Nature, 

472, 31 

White, R.E., & Shawl, S.J. 1987, ApJ, 317, 246 

Zackrisson, E., Bergvall, N., & Ostlin, C, 2005, A&A, 435, 29 

Zinn, R., 1985, ApJ, 293, 424 

Zwaan, M.A., van der Hulst, J.M., de Blok, W.J.C, & McGaugh, 
S.S. 1995, MNRAS, 273, L35 



Villcgas et al. 



APPENDIX 

NOTES ON INDIVIDUAL GALAXIES 

ugcOO '477 and ugc06138. These are the two galaxies for which we used a wide selection of control fields to estimate 
the contamination by background galaxies, due to their location at high Galactic latitude, which also contributes 
to minimize the uncertainty due to reddening conditions. The distance to the galaxies makes prohibitive a 
spectroscopic follow up with 8-m class telescopes. 

ugc03459. This is the galaxy with the most uncertain results. Its field of view is located directly through the 
galactic plane, at a latitude of b = —3.92, where any local reddening feature will produce dramatic changes in 
the results. The almost null detection of candidates in this galaxy tells us that the objects are too absorbed to 
be detected. 

ugc03587 and ugcll651. These galaxies have a high number of GC candidates detected. Nevertheless, some 
of the candidates are projected along the spiral arms of the galaxies, so some contamination due to young star 
clusters is expected. We also noticed that due to the very peculiar morphology of ugc03587, SExtractor's model 
used to subtract the light of this galaxy when detecting point sources is not as good as in the other galaxies. 
These galaxies are both close enough to be followed up with ground-based spectroscopy in 8-m class telescopes. 

ugclll31. Due to its location close to the direction of the Galactic center this galaxy is affected by strong 
extinction. It is also very difficult to obtain good control fields in its case as in this region the stellar density has 
a strong dependence on the galactic longitud, thus this is the only galaxy with just one control field. 



